Super Bloch oscillations with modulated interaction 
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We study super Bloch oscillations of ultracold atoms in a shaken lattice potential, subjected to a 
harmonically modulated mean- field interaction. Usually, any interaction leads to the decay of the 
wave packet and its super Bloch oscillation. Here, we use the phases of interaction and shaking with 
respect to the free Bloch oscillation as control parameters. We find two types of long-living cases: 
(i) suppression of the immediate broadening of the wave packet, and (ii) dynamical stability of all 
degrees of freedom. The latter relies on the rather robust symmetry argument of cyclic time [Gaul 
et a/., Phys. Rev. A 84, 053627 (2011)]. 
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Experiments with ultracold atoms in optical lattices 
are an ideal testing ground for many problems of con- 
densed-matter physics [1, 2]. In addition to good mea- 
surement access, these systems offer flexible manipula- 
tion of the system parameters, which opens the way to 
new perspectives and effects. One of the most prominent 
examples is the observation of Bloch oscillations (BOs) 
of cold atomic gases [ ] and of Bose-Einstein condensates 
(BECs) [ ] in optical lattices, subjected to an external 
force F. The semiclassical explanation for BOs is the 
following: the quasimomentum hk increases linearly with 
time, but because of the dispersion relation of a tight- 
binding model with hopping amplitude J and lattice pe- 
riod a, E(k) = 2J[1 — cos(fca)], the group velocity is a 
sinusoidal function of time. The particle does not follow 
the potential gradient but stays localized and performs 
an oscillatory motion. 

By optical means or by magnetic levitation, the lat- 
tice potential can be shaken with frequency uj. This 
results in a renormalization of the hopping amplitude, 
which can even be suppressed, the so-called dynamic 
localization [ , ]. Semiclassically, a harmonic shaking 
F(t) — AFsin(cjt) causes the quasimomentum hk = 
j dt'F(t') to oscillate rapidly and to explore /c-space re- 
gions with renormalized or even negative effective mass 
m e fi(k) oc 1/ cos(fca). Time- averaging cos(ka) leads to 
the renormalization of the hopping J e ff/J = Jo(AF/cj), 
where Jo is the Bessel function of the first kind. Thus, 
depending on the strength and the frequency of the shak- 
ing, the effective hopping J e ^ can be suppressed or even 
negative, freezing or inverting the center-of-mass motion. 

The modulation of the force around a finite mean value 
leads to a superposition of the BO with a slow oscillation 
of large real-space amplitude, similar to the one shown 
in Fig. 1. This phenomenon is known as quasi BO [ ] or 
super BO (SBO) [8]. It can be explained with a semiclas- 
sical reasoning, too: the quasimomentum performs small 
oscillations around its linearly increasing mean value. 
Then, during one Bloch cycle, it spends more time in k- 



space regions with, say, positive mass than in regions with 
negative mass, which results in a drift of the wave packet 
in the direction of the force. As the relative phase be- 
tween BO and shaking changes, the time- averaged mass 
becomes negative and the drift gets reversed, which re- 
sults in an SBO at the beating frequency. 

SBOs have been mostly studied in the linear, noninter- 
acting case [9, 10]. But ultracold bosons open more inter- 
esting possibilities. Feshbach resonances can be used to 
arbitrarily change the s- wave scattering length [ ] , even 
time-dependently [12-14]. Generically, the interaction 
leads to dephasing and decay of the wave packet. How- 
ever, the interplay of modulated interactions and BOs 
has already been investigated, and an infinite family of 




FIG. 1. (Color online) Stable interacting SBO. The initial 
wave function has Gaussian shape with width ao = 20a and 
quasimomentum p(0) = hit /2a. Force (2) with Fo = 0.2J/a, 
AF = O.6F0, l/v — 4/5, and (j)F — 0; interaction (3) with go = 
1,0 = 0. (a) Real-space density |^ n (t)| 2 from the integration 
of Eq. (1). (b) Time evolution of some magnitudes of interest: 
the interaction g(t), the inverse mass term ~ cos (pa) and the 
variation of the real-space width Aa(t) = a(t) — cr(0). The 
gray circles mark the points of symmetry used in the cyclic- 
time approach. 
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(harmonic) modulations that lead to a periodic time evo- 
lution of the wave packet has been found [15-17]. Here, 
we consider SBOs in the presence of a modulated s- wave 
scattering length. This problem is more complex because 
already the linear problem contains two frequencies and 
two phases from the BO and the shaking, respectively. In 
the remainder of this Brief Report, we tackle the problem 
with a full integration of the discrete Gross-Pitaevskii 
equation and a cyclic-time argument similar to that of 
Ref. [ ], for one particular phase of the shaking. After- 
wards we study frequencies and phases not covered by the 
cyclic-time argument by means of collective coordinates 
of a Gaussian wave packet and linear stability analysis of 
the infinite wave packet. 

Model. Our starting point is the mean-field tight- 
binding equation of motion as in Refs. [15-17], but now 
the tilt F may be time dependent: 

im n = -J(* n+ i +* n -i) +F(t)an^ n +g(t)\^ n \H n . 

(1) 

The wave function \I/ n is normalized to one, implying that 
the time-dependent interaction parameter g(t) contains 
the total particle number. We choose the force as 

F(t) = F + AFcos(ut + <j> F ). (2) 

The mean value T defines the Bloch frequency ujb = 
Foa/h. This and the frequency of the modulation uj 
are the two frequencies of interest. For concreteness, we 
choose a fixed frequency ratio uj/ujb = Ijv = 4/5 for the 
rest of this Brief Report. Thus, the super Bloch period 
Tsbo is five Bloch periods. We now search for suitable 
modulations of g(t) that allow a periodic time evolution 
of (1). A constant interaction parameter leads sooner 
or later to a decay of the SBO [ ], but harmonic modu- 
lations around zero may counteract this effect [15]. We 
choose the simplest commensurate modulation 

g(t) = gosm(u B t/v + 0). (3) 

Throughout the paper and we vary only the phase <j). 
Figure 1(a) shows the time evolution of the wave function 
for <pp = cj) = 0, a particular choice of parameters for 
which SBOs turn out to be stable. 

In general the SBO decays with time. This decay is 
most conveniently observed via the broadening of the 
wave-function momentum distribution, observed via the 
inverse participation number (IPN) of Fourier modes 

IPN = £|* fc | 4 . (4) 

k 

A decrease of this number means momentum broadening 
and decay of the wave packet. Figure 2 shows the varia- 
tion A IPN (T) = IPN (T) - IPN (0) after a certain time T 
for (j)F = 0, — 7r/2 as a function of the interaction phase 
(j). For most of the phases, the IPN goes down, meaning 
that the wave packet broadens. Figure 2(a) shows that 
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FIG. 2. (Color online) Relative change of the momentum 
IPN as a function of the phase cj) for (a) <j)F — and 
(b) 4>f = — 7r/2. Discs (with lines as guide to the eye): results 
from integration of Eq. (1) after integration times T = Tsbo, 
15Tsbo, and 30Tsbo (from top to bottom). Solid line: Result 
from the collective coordinates theory (11). Regardless of a 
possible initial increase, the IPN always goes down in the long 
run, except for those points where A IPN = 0. 

for <j)p = the phase = from Fig. 1 is one of the 
two stable points, where the IPN stays constant. Simi- 
larly, in the case = —tt/2 of Fig. 2(b), two "stable" 
points arise as well. However, their stability behavior is 
different and will be explained in what follows. 

Cyclic-time formalism for SBOs. We express the dis- 
crete wave function \I/ n of (1) in terms of a continuous 
wave function A(z,i): 

9 n (t) = A(na - x(t),t)e l ^ t)na/h+im . (5) 

With the semiclassical equations p — —F(t), x — v g = 
2 J a sin (pa/ h)/h, and hep = 2 J cos (pa/ ft), the envelope 
function obeys the nonlinear Schrodinger equation 

ihA = -Jcos(pa/K)A" + g(t)\A\ 2 A , (6) 

where higher derivatives of A have been neglected. The 
factor in front of the second derivative takes the role of 
the mass term J cos (pa/K) = h 2 /2m(t). Taking into ac- 
count the commensurability of the frequencies uo and Co>b ? 
we express the time dependence in terms of the slower 
time r = u&t/v, which increases by 2tt during one super 
Bloch period Tsbo = After integrating the force 

(2), we find the inverse mass term as 

cos(pa/h) = sin \ vr + 9l j^- [sin(Zr + (j)p) — sin^i?]} , 

(7) 

where we have chosen the initial condition such that 
ap(r = 0) = irh/2. This can always be achieved by 
choosing the origin of time. 

We know from previous works [15-17] that BOs can 
exist in the presence of suitably tuned AC interactions. 
Now, we construct a similar cyclic-time formalism for 
SBOs. The basic idea is to separate all terms of Eq. 
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(6) into the form 77/(77), where 77 is a harmonic function 
of r and / an arbitrary function. Then 77 is eliminated 
and Eq. (6) is solved as a function of 77 only. As 77 is a 
periodic function of t, the time evolution of A is periodic 
too, and we have indeed found a case of stable SBOs in 
the presence of a non-zero interaction g(t). 

The inverse mass term (7) can only be brought into 
the desired form 77/(77) if = [or any phase shift that 
allows for a common node of sin^r and sin Zr]. Then, 



cos(pa/H) = sin(r)J r (cosr), 



(8) 



which is the most general 27r-periodic function that re- 
spects the odd symmetry with respect to r = mod n. 
This determines 77 = cost. Thus, the left hand side of 
(6) is A = fjd^A. Finally, we may choose any modu- 
lated interaction of the form g = 77*7(77). In particular, 
g(t) = g s'm(F t/v), i.e., Eq. (3) with = fulfills the 
cyclic-time condition. The symmetries involved in this 
argument can be observed in the lower panel of Fig. 1: 
g(t) and cos(p(t)a/H) share the nodes of 77, marked with 
gray circles, and are odd with respect to these points. 
All physical quantities, e.g. the width variation Act, are 
functions of 77 only and are thus even with respect to the 
mentioned points. 

So far, we have explained the stability observed in 
Fig. 1 and the stable points of Fig. 2(a). It turns out 
that the phase = —tt/2 of shaking considered in 
Fig. 2(b) cannot be written in the form required by the 
above cyclic-time argument, because the inverse mass 
term cos ap/H does not exhibit points with odd symme- 
try. Nevertheless, there are two particular phase shifts, 
(j> ~ — 0.267T and (j) ~ 0.747T, where the broadening of the 
wave packet is suppressed. One can understand this us- 
ing a collective- coordinates theory. The idea is to reduce 
the complexity of the problem by parametrizing the wave 
function by only two coordinates, its position x and its 
width y/w [15, 18]: 



*»(') = ^ 



w 



ipna/h-\-ib(na—x) /h 



(9) 



A(u) is a continuous normalized Gaussian wave function 
with variance one. The equations of motion for the collec- 
tive coordinates, x and w, and their conjugate momenta, 
p and 6, are 
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h x 

J a 
J 

hw 



b = 



= ~F(t) , 
= 2a sin (pa/ K) 
I6w 2 b 2 /ti 



1/4 + 4b 2 w 2 /H 



2w/a 2 



4w 2 /a 4 



cos(pa/h) 



8^Jw 3 / 2 ' 



wb 



8 — cos(pa/h) . 



(10a) 
(10b) 

(10c) 
(lOd) 



J a 2 ~ h 

A necessary condition for stability within the collective- 
coordinates approximation is that w(t) returns to its ini- 
tial value after a full super Bloch period Tsbo- In this 
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FIG. 3. (Color online) Momentum- space portraits in the ref- 
erence frame of the center of momentum p for shaking (2) and 
interaction (3), with a = 20a, F = 0.2 J/a, AT = 0.6T , 
go = 1. (a) Unstable sine shaking (j)F — —tt/2 with = 0; 
(b) sine shaking with phase = — 0.267T adjusted according 
to Eq. (11); (c) stable case of Fig. 1, (\>f — 0, (j) — 0. 



regard, we can compute analytical results in the limit 
of a wide wave packet w(t) ~ <Jq ^> 1. Then, Eqs. (3) 
and (10c) give b(t) ~ — avg /[8^o-qUJb] cos(o;b^/^ + 4 > )- 



With (lOd), and assuming that 



1, this yields 



dt 



w 



{l-v) 



Tsbo 



y/7ra 3 Jgo 

o-l(hw B /yf 



sin ( 



(ii) 

with (j)Q = (I — v) \(j)F + aAi ? sin(0i?)//kj]. Depending on 
the phase (j) of the interaction, the wave packet contracts 
or spreads after a full super Bloch cycle Tsbo- This is 
also reflected in the IPN shown in Fig. 2. The dynamics 
can be strictly periodic only for = O mod tt. With 
the values I = 4, v = 5, <j)p = —tt/2, and aAF = 0.75fkj, 
this yields </>o = 0.747T, in agreement with Fig. 2(b). 

In Fig. 3, the momentum-space portrait of SBOs is 
shown for different phases of the shaking and of the inter- 
action. Panel (a) shows an unstable case of sine shaking 
with immediate broadening of the momentum distribu- 
tion. In panel (b) the phase of the interaction is adjusted 
as required by the collective-coordinates criterion of Eq. 
(11) and it shows the behavior presented by the two ap- 
parently "stable" points in panel (b) of Fig. 2. Here, 
the immediate broadening is indeed suppressed; but on a 
longer time scale (~ 20Tsbo = 100Tb), side peaks grow 
far from the central wave packet. This is an indicator for 
a dynamical instability and does not happen in the stable 
case with cosine shaking (see Fig. 1), whose /c-space por- 
trait is shown in Fig. 3(c) for comparison. In this case, 
there is neither immediate broadening nor an instability. 
The extremely small residual broadening, is due to effects 
beyond the approximation made in Eq. (6). 

Linear stability analysis. Collective coordinates rely 
on a smooth Gaussian profile of the wave function and 
cannot describe features on short length scales. We close 
this gap with a linear stability analysis of the infinitely 
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FIG. 4. (Color online) Stability prediction A q (black: stable, 
white: unstable) for fluctuation mode q = k — p as a function 
of the interaction phase <j> for cosine shaking <j)p = (a) and 
sine shaking (\)f — —tv/2 (b). The dashed green line indicates 
the phase (j) g taken in Figs. 3 (b) and (c), respectively. 



extended Bloch oscillating wave packet. The method is 
the same as in Ref. [15], only that cos(Ft) is replaced with 
cos(pa/h) from Eq. (7) and the period with the super 
Bloch period. For each plane-wave modulation q of the 
wave function, the monodromy matrix is obtained by in- 
tegration of the (linearized) equations of motion over one 
period with two different initial conditions. The mode q 
is stable if A q , which is half the trace of the monodromy 
matrix, is less or equal to one. Conversely, it is unstable 
if \A q \ > 1 [ ]. Figure 4 shows maps of stability for both 
cases of primary interest, cosine shaking ((pp = 0) and 
sine shaking (<pp = — 7r/2). In the first case, for values 
4> = rar, n G Z, all modes q remain stable, in concor- 
dance with the cyclic-time argument and with Fig. 2(a). 
In the second case, the stable phase of collective coordi- 
nates agrees with the stability for small values of q only, 
and for larger values instabilities occur. Notably, the first 
region of instability is located around q = 0.5. This is 
just the region where the instability occurs in the full 
integration, as shown in Fig. 3(b). 

Conclusions. We have studied the stability of com- 
mensurate super Bloch oscillations under a time- 
dependent interaction parameter that is modulated with 
the super Bloch frequency. The phase of the interaction 
can always be adjusted such that the direct broadening 
of the width is suppressed [see Fig. 2 and Eq. (11)]. But 
there are more degrees of freedom, which may be subject 
to dynamical instability [Fig. 3(b)]. These can be stabi- 
lized by fulfilling the cyclic-time argument [Fig. 3(c)]. 

The dynamics of BOs and SBOs strongly depends on 
the relative phase evolution between neighboring sites 
[4, 8]. This means that interactions distort the density 
profile of the driven Bose-Einstein condensate and lead 



to the destruction of such oscillations in general. There- 
fore, our proposal to tune the atomic interaction time 
dependently is important to reduce this effect and to al- 
low for the detection of such oscillatory dynamics of a 
condensate. The results of this work can be considered 
to improve the quality of SBOs in order to be used to en- 
gineer matter-wave transport over macroscopic distances 
in lattice potentials in a more reliable way, which should 
also be relevant for atom interferometry [20]. 
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